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ABSTRACT 

Astronomical seeing is quantified by a single parameter, turbulence integral, in 
the framework of the Kolmogorov turbulence model. This parameter can be routinely 
measured by a Differential Image Motion Monitor, DIMM. A new instrument, Multi- 
Aperture Scintillation Sensor (MASS), permits to measure the seeing in the free atmo- 
sphere above ~0.5km and, together with a DIMM, to estimate the ground-layer seeing. 
The absolute accuracy of both methods is studied here using analytical theory, numerical 
simulation, and experiments. A modification of the MASS data processing to compensate 
for partially saturated scintillation is developed. We find that the DIMM can be severely 
biased by optical aberrations (e.g. defocus) and propagation. Seeing measurements with 
DIMM and MASS can reach absolute accuracy of 10% when their biases are carefully 
controlled. Pushing this limit to 1 % appears unrealistic because the seeing itself is just a 
model-dependent parameter of a non-stationary random process. 
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1 INTRODUCTION 

Measurements of astronomical "seeing" are performed for se- 
lecting new sites and supporting operation of existing tele- 
scopes. Recently, a st andard Differential Image Motion Mon- 
itor (DIMM) method jMartid 1 19871 : ISarazin & RoddieJ 1 19901) 
has been complemented with a new technique, Multi-Aperture 
Scintillation Sensor (MASS) (Komilov et al. 2003). This new 
instrument is based on the analysis of scintillation and permits 
to measure the seeing in the free atmosphere, isoplanatic angle, 
and adaptive-optics (AO) time constant. Both MASS and DIMM 
require only a small telescope and can be combined in a single 
instrument ( Kornilov et al. 2007). MASS-DIMM site monitors 
gradually become a new standard. 

We can evaluate the ground-layer (GL) seeing produced in 
the first 0.5 km above the observatory by subtracting the turbu- 
lence integrals measured with DIMM and MASS. It is impor- 
tant to measure the GL seeing for evaluating the performance 
of Ground-Layer Adaptive Optics (GLAO). However, subtrac- 
tion only works when both methods are accurate, i.e. deliver 
results on the absolute scale. Absolute accuracy of seeing mea- 
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surements becomes critical for the site comparison, where the 
differences are often below 10%. 

In principle, both DIMM and MASS should give accurate 
results when their instrument parameters are set correctly and 
their intrinsic biases are understood and removed. Here we in- 
vestigate these biases in detail, quantify them, and propose cor- 
rections. Biases of the DIMM method have been alread y ad- 
dressed in the literature l Martinlll987l : iTokovinid 12002 j) . We 
continue by considering additional effects such as light propaga- 
tion and optical aberrations. The analysis of the MASS method 
given by Tokovinin et alj l l2003l) is extended by studying small 
departures from the weak-scintillation theory which cause a sys- 
tematic bias, "over-shoots", and can be corrected by a modified 
data processing. 

This work has been stimulated by the need to get accurate 
results from the existing suite of MASS-DIMM instruments, de- 
scribed in the accompanying paper (Kornilov et al. 2007). Cor- 
rect setting and operation of these instruments, also critical for 
getting accurate data, is addressed in that paper. Here we con- 
centrate on the theoretical analysis and simulations, trying to 
formulate practical recommendations and recipes while avoid- 
ing mathematical complexity as much as possible. 

We begin by asking the question "what is seeing and how 
accurately can it be defined and measured?" in Sect.|2] In Sect.[3] 
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Table 1. Quantities relevant to seeing 



Plane wave 



Phase screen 



Quantity 


Units 


Seeing (FWHM) 


rad 


Turbulence integral 


ml/3 


Fried parameter 


m 


Phase power spectrum 


m2 



Formula 



eo = 0.98A/ro 

ro = [0.423(27r/A)2 J]-3/5 
ro = l.OlA/eo 

*^{/) = ci/r"/^ 

C = 0.00969{27r/A)2 J 



C = 0.0229r, 



-5/3 







the MASS method is studied under conditions of reaUstic (not 
weak) scintillation. Then in Sect. |4] we address two previously 
neglected effects in a DIMM - propagation and optical aberra- 
tions - and show that they can cause a significant bias, in ad- 
dition to the known DIMM biases. Our conclusions and recom- 
mendations are formulated in Sect.|5l 



2 WHAT IS MEASURED BY A SEEING MONITOR? 
2.1 Turbulence parameters 

Turbulence measurements are based on the standard theory 
jTatarskiilll96ll : lRoddieilll98lh . For convenience, major atmo- 
spheric parameters are recalled in Table [T] These definitions as- 
sume a scale-free Kolmogorov turbulence spectrum with only 
one parameter, turbulence strength. This single parameter can be 
expressed equivalently by eo, '"o, or J. These quantities (except 
J) depend on the imaging wavelength A which is assumed here 
to be 500 nm if not specified. Accurate measurement of this sin- 
gle parameter (called "seeing" in a broad sense) is the purpose of 
a seeing monitor. The atmosphere can be split into an arbitrary 
number of zones (or layers), and we also want to measure the 
seeing produced by each of the layers - the turbulence profile. 
Other atmospheric parameters (time constant, isoplanatic angle), 
not considered here, are also of interest to modern astronomy. 

The seeing eo is often considered to be the angular image 
spread caused by the atmosphere. This interpretation is only ap- 
proximate because even in a perfect telescope the actual long- 
exposure point spread function (PSF) depends on several addi- 
tional parameters (outer scale, wavelength, telescope diameter, 
guiding). Current seeing monitors measure only one number eo 
which is necessary, but not sufficient for accurate prediction of 
the atmospheric PSF or other turbulence-related optical quanti- 
ties. 

It is always implicitly assumed that the statistical properties 
of turbulence are stationary, while in fact they are not. Any mea- 
surement refers only to the particular moment in time and to the 
particular viewing direction. The non-stationarity precludes very 
precise measurements of atmospheric parameters because aver- 
aging over a larger spatial or temporal sample of turbulence does 
not lead to the improved statistical precision. This aspect, often 
overlooked, also makes seeing measurements intrinsically irre- 
producible. Comparisons between seeing monitors or between 
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Figure 1. Schematic representation of a seeing monitor. The exclamation 
marks and text in italics show approximations or uncertainties present in 
any seeing measurements. 



a monitor and a telescope can be only made in a statistical 
sense, with associated non-stationarity errors. Fortunately, in a 
MASS-DIMM instrument both channels sample the same turbu- 
lent path, hence non-stationarity does not affect the precision of 
the GL seeing estimate. 

The seeing is not a well-defined physical quantity like 
length or mass, so the intrinsic accuracy of seeing monitors can- 
not be arbitrarily high. It is unrealistic to expect a relative ac- 
curacy better than 1% because the "seeing" cannot be defined 
with such a high accuracy. It is shown below that keeping biases 
within few percent is not easy. 



2.2 Seeing monitors 

A seeing monitor measures some statistical properties of phase 
and/or amplitude of a light wave at the ground and interprets 
them in terms of the model parameter, seeing. Several approxi- 
mations are involved in this process, as illustrated in Fig.[T]and 
discussed below. 

Both DIMM and MASS are sensitive to phase distortions 
with spatial scales below 1 m where the Kolmogorov model 
works well. A DIMM with small apertures is not affected by 
the finite turbulence o uter scale Lp (less than 1% bias on vari- 
ance for Lo > 4 m, cf iBorgnino, Martin. & Ziaj[T993) . On the 
other hand, the absolute image motion in a small 10-cm tele- 
scope is influenced by a finite Lo, typically at ~ 10% level. Thus, 
a site-testing instrument based on the absolute (non-differential) 
image motion gives a biased ro estimate if Lo is not known. In 
practise, such seeing monitors are no longer used because they 
are also affected by mechanical errors (wind shake, tracking). 
For the same reasons, the theore tically perfect interferom etric 
method of seeing measurements ( iDaintv & Scaddanlll975l) has 
not become widely adopted. 
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Figure 2. Intensity screens of 1 size produced by a single turbulent 
layer and 10 km propagation in conditions of weak scintillation (left, 0'.'2 
seeing) and strong scintillation (right, 1" seeing). Wavelength 0.45 fim. 

2.3 Propagation 

Seeing measurements are affected by the light propagation trans- 
forming pure phase distortions to a mixture of phase and ampli- 
tude distributions. The spatial spectra of the light phase (p and 
logarithm of amplitude x after passing through a weak turbulent 
layer and propagation over a distance z are 

$^(/) = 0.0229 r-"/'j/r"/^cos2(^A2|/l"), (1) 
= 0.0229 ro~'/'|/r"^'sin^Mz|/|2). (2) 

The spectrum of intensity fluctuations (scintillation) is <I>7 — 
4:^x- The amplitude and phase are not correlated at any given 
point, but their cross-spectrum is not zero at / 7^ 0, being pro- 
portional to the product of the sine and cosine terms. 
The scintillation index is defined as 

s= = {AI^)/{lf, (3) 

where I is the instantaneous light intensity received through 
some aperture, AI is its fluctuation. In the small-perturbation 
regime, s'^ <^ 1, formula ([3} is equivalent to the variance of 
the log/ = 2x- The scintillation index is calculated then by 
integrating the amplitude spectrum I2J with a suitable aperture 
filter. The effect of several turbulent layers is simply additive. 

The equations ([TJ and (|2} are only approximate for a real 
(not weak) turbulence. Even near the zenith, scintillation indices 
> 1 were measured. This regime of strong scintillation cor- 
responds to the onset of focusing when the deviations from the 
standard theory become quite significant. Figures [2] and [3] show 
simulated scintillation signals and their spatial power spectra for 
weak and strong scintillation (cf. Sect. I2.4t . In the latter case, 
the pattern is dominated by small spikes due to focusing of the 
light by the "lenses" created by the high-atmosphere turbulence. 
The size of these spikes is less than the size of the lenses and, ac- 
cordingly, the power at high spatial frequencies increases com- 
pared to the weak-scintillation spectrum while at low fre- 
quencies it decreases. The appearance of extra power can be ex- 
plained qualitatively in the following way. Under weak scintilla- 
tion, each sinusoidal component of the phase screen creates in- 
tensity modulation with the same frequency. As the amplitude of 
the phase perturbation increases, the phase gratings start to pro- 
duce second and higher-order or crossed harmonics in the inten- 
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Figure 3. Radially-averaged power spectra of scintillation correspond- 
ing to the two cases of Fig.|2] "Noisy" curves - simulated, smooth curves 
- weak-perturbation theory, Eq.|2] 

sity distribution at the ground. Individual layers are no longer in- 
dependent but cross-modulate each other. The phase is affected 
by the saturation in a similar way as the amplitude. 

2.4 Simulation tool 

Given the lack of accurate wave-propagation theory, the best 
way to study higher-order propagation effects is by numerical 
simulation. Simulations are also indispensable in evaluating var- 
ious instrumental effects. 

We use the Fourier method of generating random phase 
screens with Kolmogorov statistics. A 2-dimensional array of 
zero-mean Gaussian complex random numbers is created, their 
amplitudes oc (/ is the spatial frequency, C is the 

coefficient from Table [TJ and phases distributed uniformly in 
the interval (— vr, n). The Fourier transform (FT) of this array 
creates the phase screen. It is well known that this method un- 
derestimates the low-frequency components of turbulence, and, 
notably, produces wrong phase structure functions. The reason 
is that any function obtained by discrete FT is periodic, with a 
period equal to the grid size. However, the effects of the numer- 
ical outer scale can be neglected if the aperture size is a small 
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fraction of the screen size. Indeed, we checked that in the weak- 
perturbation regime our numerical simulations reproduce the an- 
alytical scintillation spectra and DIMM response to within 2%. 

The propagation of wave-fronts is calculated by the spec- 
tral technique. If Ui{r) is the amplitude of the light wave before 
propagation, and Ui{f) is its FT, then the FT of the amplitude 
U2{f) after propagation over a distance z is obtained by fre- 
quency filtering: 

f>2(/) = (7i(/)exp(-i^^Aj/|2). (4) 

This method is computationally fast, involving only two 
FTs. Its drawback is that in fact it describes the propagation 
in a waveguide with a cross-section equal to the simulation 
grid and reflective walls. However, the propagation of periodic 
phase screens does not present problems near the grid bound- 
aries because the amplitude "wraps around" the simulated do- 
main boundaries and does not produce artifacts. 

We generate large screens of complex light amplitudes at 
the ground resulting from the propagation through one or sev- 
eral turbulent layers at various altitudes. Typically, with 1024^ 
points and 0.5 cm sampling the grid size is 5.12 m. To simulate 
the data sequence of a seeing monitor, the screens are shifted in 
both coordinates by Vt, where V is the wind speed and r is 
the sampling time interval. The shift is directed at some small 
angle with respect to the x-axis. In this way, the aperture moves 
mostly in x but is displaced in y by ~ 0.2 m per line, eventually 
covering the whole screen one or several times. There are no 
adverse effects at the borders because the screens are periodic. 
The statistical averaging is sufficient to simulate typical 1-min 
measurements. 

The complex amplitude of light at the aperture of a seeing 
monitor is re-sampled on a finer grid and used to calculate the 
measured quantities such as the spot images in DIMM or fluxes 
in MASS. The effects of finite CCD pixels and detector noise 
can be studied as well. The current simulator has some limita- 
tions. The wind speed is common to all layers. The blur during 
a finite exposure time in DIMM is not simulated. In MASS, we 
simulate the exposure time r by a linear blur of the apertures 
over a distance Vt. The simulation is usually monochromatic. 



3 ACCURACY OF THE MASS METHOD 
3.1 From scintillation to seeing 

The MASS instrument is based on the spatial analysis of inten- 
sity fluctuations at the ground level. The spatial scale of the scin- 
tillation "speckle" produced by turbulence at a distance z is of 
the order of the Fresn el radius rp = ^/Xz, i.e. ~ 10 cm for a 
10-km propagation ( Roddiedll98l l). The spectrum l|2j reaches 
maximum at the spatial frequency |/| ~ r-p ^. 

Light from a bright star is detected in MASS with four con- 
centric annular apertures with diameters from 2 cm (inner) to 
8 cm (outer). The size of these apertures is of the order of rp 
and they act jointly as a spatial filter, permitting to dis-entangle 
scintillation originating at different altitudes. This is achieved in 
several steps. 



Step 1. The sequences of photon counts from 4 apertures 
with individual micro-exposure of 1 ms are acquired and pro- 
cessed to calculate 10 scintillation indices - normal indices s\ 
for each aperture A and 6 differential indices s^b for pairs of 
apertures A and B. The formulae for calculating the indices and 
subtracting the bias c aused by the photon noise are given in 
l lTokovinin et al.ll2003l) . The indices s\ and s^b are equal to 
the variance of the natural logarithms log Ja and log(JA//B) in 
the limit of small fluctuations of the light intensities I a and 7b, 
<^ 1. The turbulence theory usually operates with the vari- 
ance of the logarithm. However, we measure the photon counts 
which can be zero and fluctuate even at constant light, hence the 
indices should be calculated from the normal (non-logarithmic) 
variances. 

Step 2. A linear relation between the observables (indices) 
and the turbulence profile C^{z) is established by the weak- 
perturbation theory, 

sl = J Wk{z) Cl{z)dz, (5) 

where the weighting function (WF) Wk{z) describes the alti- 
tude response of a given aperture or aperture combination k. 
For a weak Kolmogorov turbulence, the WF depends only on 
the aperture geometry and spectrum of detected radiation and 
can be derived from ^ (Tokovinin 2002b, 2003). A normal 
index for an aperture of diameter Da corresponds to the low- 
pass filtering of scintillation, passing j/j < 1/Da, whereas a 
differential index corresponds to a band-pass spatial filter with 
1/-Db < l/i < 1/-Da. Thus, MASS with its centimetric aper- 
tures is sensitive to the turbulence of centimetric scales where 
the Kolmogorov model is adequate. 

It has been shown that the differential index in two con- 
centric annular apertures is almost independent of the propa- 
gation distance z for z > zab ~ D\^/\, wher e -Dab = 
(-Da + _Db)/2 is the average aperture diameter UTokovininl 
l2002bl l2003[) . It means that the scintillation index gives a di- 
rect measure of the turbulence integral, hence seeing, produced 
at distances beyond zab- On the other hand, normal scintilla- 
tion indices increase as z^ , with /3 = 5/6 for smal l distance 
2 < -D^/A and /? = 2 for z > -D^/A ( lRoddieJll98ll) . Accord- 
ingly, most of the scintillation is produced by the high layers. 

Step 3. Using the known WFs, the set of 10 indices is fitted 
with a model of 6 thin turbulent layers at altitudes hi of 0.5, 1, 
2, 4, 8, and 16 km, with turbulenc e integrals Ji in each layer as 
parameters jTokovinin et al.ll2003h . The zenith angle 7 is taken 
into account by the model, z = h sec 7. In reality, the turbu- 
lence is distributed continuously in altitude with a profile (h) . 
The integrals Ji delivered by MASS are approximately equal to 
J Cn{h)Ri{h)dh, where the response functio ns Rj (h) resem- 
ble tr iangles in log h coordinate centred on hi jTokovinin et al.l 
l2003h . The sum of all Ri (h) is close to one for h > 0.5 km. 



3.2 Over-shoots and their correction 

MASS relies on the small-perturbation theory, assuming that the 
scintillation is weak, s'^ ^ 1, and that the combined effect of 



© 2007 RAS, MNRAS OOO.fTlfn] 



Accurate seeing measurements with MASS and DIMM 5 



1 .4 - 
< 1.2 - 

QJ 

~o 

o 

^ 1 .0 - 

I 

0.8 - 

0.6 - 
0.01 

1 .4 - 




m 1 9 



2 1.0-S[ffl- 



0.8 - 



0.6 



0.01 



1.4 



1 .0 



- H 



- ^ 0.8 



0.10 

Scintillation index A 



1 .00 



0.6 - 

0.01 

1 .4 



O 



1 .2 



1 .0 



0.8 



0.6 



0.10 

Scintillation index A 



1.00 0.01 



0.10 

Scintillation index A 



1 .00 



'^7 



0.10 

Scintillation index A 



1 .00 



Figure 4. Ratio of measured to modelled scintillation indices. Selected normal (top row) and differential (bottom row) indices in some MASS apertures 
as determined from the simulations without coiTection (s^, crosses) or after connecting by Eq.|6](s^, empty squares) are divided by the weak-scintillation 
indices s^. The dotted vertical line marks the correction limit = 0.7. 



several turbulent layers is additive (Eq.|5}. In practise, scintilla- 
tion is not always weak (Sect. [23l l. In this case, profile restora- 
tion by the linear method, as implemented in MASS, leads to the 
over-estimated total turbulence integral (free-atmosphere see- 
ing) and to the shift of the restored profile to lower altitudes. 

Although the theory of strong scintillation ha s been ad- 
dressed in many papers (e.g. [Andrews et al.l 1 19991 and refer- 
ences therein), there is no quantitative description of the inten- 
sity power spectrum available. In order to extend MASS oper- 
ation to moderately strong scintillation, we rely exclusively on 
numerical simulations (Sect. l2.4l (. Poly-chromatic light was sim- 
ulated by an equal mix of three wavelengths of 0.4, 0.45, and 
0.55 fim. One or two phase screens at various altitudes were 
simulated. The physical size of simulated phase screens is 2.5 m 
across, with 2.5 mm pixels. The power spectrum of intensity 
fluctuations at the ground was computed (e.g. Fig. [3} and con- 
verted into normal and differential scintillation indices of MASS 
by integrating it with suitable spatial filters. 

The values of the measured indices inform us on the 
strength of the scintillation. Our approach is to correct the mea- 
sured indices semi-empirically and to bring them into agreement 
with the weak-scintillation theory. Let be the 10-element vec- 



tor of measured indices, Sq - the vector of theoretical indices 
expected in the linear theory without saturation, and - the 
vector of corrected indices. A rather general correction formula 
can be written as 



The rationale for selecting this formula is that it automatically 
removes correction for weak scintillation, and that this type of 
formula works well for the differential indices. Here, Z is the 
10 X 10 correction matrix. 

We determine the correction matrix Z empirically from the 
results of simulations, by least-squares fitting. The fitting is re- 
stricted to the relevant scintillation range 0.1 < si < 0.7 be- 
cause weak scintillation does not need correction (but adds sta- 
tistical noise), while correcting stronger scintillation is hopeless. 
The matrix Z is found by least squares as 

Z= (SS^)-'(S^Y), Y = S/So-l, (7) 

where the matrices S, So, Y have dimensions 10 x M, M being 
the number of simulated cases. We simulated both single and 
double layers at various altitudes with various strengths of each 
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layer. All results were mixed together in calculating the matrix 
Z, with a total of M — 35 cases satisfying the condition 0.1 < 
si < 0.7. 

Figure |4] demonstrates the success of this approach with 
plots of the simulated (pluses) and corrected (squares) indices 
divided by the weak-scintillation ones. The average ratio sJ/sq 
of all indices differs from 1 by less than 1% after correction, 
the scatter is also significantly reduced. The best correction is 
achieved for s\^, rms s'^/sq scatter 1.6%, and the worst correc- 
tion - for Sq, rms scatter 4.6%. 

The correction of indices with the matrix Z "learned" from 
the simulations is implemented in the current version of the 
MASS software, TURBINA. Old data can be re-processed with 
this program. An example of successive "overshoot" correction 
is shown in Fig.|5] Our empirical correction technique works for 
s% < 0.7. The correction matrix is determined only for the typ- 
ical aperture diameters and bandpass used in MASS, it has to be 
revised if these parameters change. 

A simplified method of correcting MASS seeing for over- 
shoots has been established earlier and works quite well. If 
'''o,MASS is the Fried parameter of the free-atmosphere seeing 
measured with the old MASS software, then the corrected ro 
will be 



ro 



To, MASS (1 + 0.7sa)°'' 



(8) 



3.3 Temporal sampling 

The photon counts in MASS are sampled with the exposure time 
At = 1 ms. Averaging the signals during At reduces the fluctu- 
ations, causing a bias in the measured indices. This bias is cor- 
rected in the software by calculating the indices with exposure 
time At and 2At and extrapolating linearly to zero exposure. 
This extrapolation is calculated as Sg = 1.5si — 0.5pi, where 
Si is the measured ind ex and pi is the covari ance with a time lag 
of 1 sampling period (Toko vinin et al.ll2003h . 

In October 2004 we recorded MASS signals with a faster 
sampling At = 0.25 ms under rapid- turbulence conditions. It 
turned out that the linear extrapolation actually over-corrects the 
indices s\ a nd Sar by as m uch as ~ 6%. An analytical study 
made earlier l lTokovinirj2002 b) also concluded that a less drastic 
correction works better. Accordingly, in the new MASS software 
the correction is halved, Sq = 1.25si — 0.25pi. The remaining 
bias is under 2% even on the fastest indices s\ and s^g. 



3.4 Inner turbulence scale 

Turbulence spectrum at high spatial frequencies (comparable to 
the turbulence inner scale) may have excess power compared to 
the K olmogorov model, the so-called Hill bump tAndrews et al. 
Il999h . The inner scale is usually of the order of few millimetres, 
but it can reach centimetric values in the upper atmosphere. This 
phenomenon can potentially increase the scintillation and lead 
to over-shoots. 

Apart from the extra scintillation, small distortions unac- 
counted for by the Kolmogorov model must cause excessive 
image blur, appearing as an extra halo of the PSF. We tried to 
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Figure 5. Example of MASS data re-processing. Turbulence profiles at 
Cerro Tololo for the night of September 30, 2004 are plotted as bars 
with length proportional to the integrals Ji (cf. the scale on the right). 
Top - original results, middle - with over-shoot con'ection. The free- 
atmosphere seeing calculated from these profiles is compared on the 
lower plot. 
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Figure 6. Comparison of the optical transfer function of DIMM spots 
(full line) with the models based on the Kolmogorov spectrum and the 
seeing measured by the DIMM (dashed line) and MASS (dash-dotted 
line) channels. Cerro Tololo, October I, 2004 2:58UT. The dotted hne is 
a diffraction-limited OTF. 

detect such effect by analysing average re-centred images of 
a star in the DIMM channel of a MASS-DIMM. With 10-cm 
DIMM apertures, the spots are always diffraction-limited under 
good seeing, but become broadened as the seeing degrades. Fig- 
ure |6] shows an example of the 1 -dimensional optical transfer 
function (OTF) of average re-centred DIMM spots in spatial co- 
ordinate r = A/. It is compared to the product of the short- 
exposure atmospheric and diffraction-limited TFs calculated for 
A — 0.55 iim with ro measured by both MASS and DIMM. 
We do not notice any departure from the DIMM model at cen- 
timetric scales and conclude that the 1.18" seeing adequately 
described the spot profile, without any perceptible effects of the 
"Hill bump". On the other hand, MASS over-estimated the see- 
ing because of the un-corrected over-shoots. This analysis was 
repeated several times on other nights. We conclude that the ef- 
fect of the finite inner scale on MASS can be safely neglected. 



propagation through turbulence was calculated with I cm sam- 
pling (Sect. I2.4t . The same amplitude screens can be re-used 
with varying DIMM parameters (e.g. aberrations). The simu- 
lated baseline is 0.24 m (even number of pixels). The centroid 
window follows the spot, as in typical DIMM instruments, be- 
cause otherwise the spots would move with respect to the win- 
dow and the response under bad seeing would be diminished. 
With 4000 spot pairs in a typical simulation, the statistical error 
of the differential variance (hence of the derived response coeffi- 
cients) is 1.6%. Apart from the statistical errors, subtle details of 
the algorithm and its implementation may influence the results 
of simulations, which are not exact in the absolute sense. 

4.1 Response coefficient of an ideal DIMM 

In a DIMM, two circular portions of the wavefront are isolated. 
The variance of the differential wave-front tilts in longitudinal 
(parallel to the base) a'f and t ransverse directions is related 
to the Fried parameter ro as l lMartin|[l987l :ls arazin & Roddieil 
[l990; Tokovinin 2002i) 

al, = Ki,t {\/Df {D/rof\ (9) 

The response coefficients of DIMM Ki and Kt depend on 
the B / D ratio and on the kind of the tilt measured. Usually, the 
tilt is evaluated from the centroids of two stellar images formed 
by the sub-apertures; in this case it corresponds to the wave- 
front gradient (G-tilt). The response coefficients for the G-tilt as 
a function of 6 = B/ D can be approximated by the formulae 
lTokoviniij|2002al) 

Ki = 0.340 (1 - 0.5706"^''^ - QMQh-'^'^) 

ift = 0.340 (1 - 0.8556"'/^ +0.0306^^/^). ^ ' 

For our example, h = 2.5, Ki = 0.1956, and Kt = 0.1270. 

The ro parameter (and seeing) is computed from the mea- 
sured differential image-motion variance by inverting Eq.|9] 
Hence, it is proportional to the response coefficients to the power 
3/5. A 10% error in the response entails a 6% error in seeing. 



4 ACCURACY OF THE DIMM METHOD 

DIMM has a reputation of a "fool-proof" technique for mea- 
suring seeing. However, this is not true because its results are 
affected by optical propagation, centroid algorithm, aberrations, 
exposure time, noise, etc. Some of these effects can cause a bias 
much larger than 10%. The accuracy of the DIMM method is 
studied below by both numerical simulations and analytical the- 
ory. The analytics is useful for understanding the small-signal 
response of a DIMM which turns out to be rather different from 
the usual assumption that DIMM measures tilts. On the other 
hand, simulations reveal the faults of the weak-perturbation the- 
ory, as for the MASS. 

We consider here a DIMM instrument with typical parame- 
ters: aperture diameter D = 10 cm and baseline B = 25 cm. In 
the numerical simulations, we selected a pixel scale of 0'.'32 and 
a field size of 10" around each spot. The complex amplitude of 
the monochromatic (A = 0.5 fim) light waves at the ground after 



4.2 Response of the centroid algorithm 

In a real DIMM instrument, the image motion is estimated by 
calculating centroids of the spots. Only a sub-set of detector pix- 
els is used in order to reduce the influence of the noise. These 
pixels are selected either by setting a threshold well above the 
background noise or by defining a window around the brightest 
pixel. Both approaches can be expressed by a general formula 

Cx ^''^WijXijIij/Io, /o = ^ (11) 

where Cx is the estimated centroid a;-coordinate, lij are pixel 
intensities, Xij are their a;-coordinates. The weights Wij equal 
one for selected pixels and zero otherwise, although a more so- 
phisticated weighting scheme could be adopted. Here we explore 
circular windows where w — 1 for pixels at a distance less than 
S\/D from the spot centre. The parameter S defines the radius 
of the centroid window in units of the diffraction spot size X/D. 
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0.0 

X — coordinate 

Figure 7. Cuts through the centroid filter functions F^{x) in the x- 
direction are plotted in ai'bitrary units against normalised ^'-coordinate, 
where ±1 corresponds to the pupil border. Two centroid windows: nar- 
row 6 = 1 and wide (5 = 5. No aberrations, D = 0.1 m. 
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Figure 8. Centroid filter functions F^pix) and F^{x) for a defocus of 
1 radian rms and 5 = 5. Compare with Fig.|7] 



m ^ ^ ^ 



Formula Jilt is only an approximation to the true image 
centroid. As a result, the response coefficient of a DIMM is not 
exactly equal to its theoretical value for G-tilt, but rather depends 
on the details of centroid calculation. We developed an analyti- 
cal formula relating the distribution of the light phase and ampli- 
tude at the pupil to the change of the image centroid /S.Cx (Ap- 
pendix 0, valid for small perturbations. It is analogous to the 
Taylor expansion of PSF developed for high-contrast imaging 
jPerrin elTal .l2003b . The centroid shift Ac^, is equal to the sum of 
the integrals of atmospheric perturbations ^{x) and xi.^) mul- 
tiplied by \he filter functions F^{x) and F-^{x), respectively. 
Knowing the spectra of the perturbations, we calculate the small- 
signal response coefficients of a DIMM with Eq. IA9I The rela- 
tive accuracy of our calculation is 3% or less. 

Figure|7]shows the filter functions of the centroid estimator 
jilt in the pupil plane calculated by Eqs. IA7IA8I For a narrow 
(5 = 1) centroid window, the filter resembles remotely a Zernike 
tilt, while for a wide {5 = 5) window it is closer to a gradi- 
ent averaged over the aperture (opposite spikes at the edges). 
The DIMM optics is assumed perfect, therefore the centroid is 
insensitive to amplitude fluctuations at the pupil (scintillation), 
F^{x) = 0. 

The situation becomes more complicated when we consider 
a realistic DIMM instrument with some optical aberrations. In 
this case, the centroid Cx is affected by both phase and amplitude 
fluctuations at the pupil (Fig. [8j. This effect can be understood 
qualitatively: a defocused spot resembles a pupil image, hence 
the centroid estimator becomes sensitive to the gradient of pupil 
illumination. Sequences of simulated spot images (Fig.|9]l illus- 
trate this situation. In a DIMM with perfect optics, the spots are 
sharp, while their intensities fluctuate because of the scintilla- 
tion. The defocused spots look more distorted by the same see- 
ing and their centroids move more, biasing the measurements. 

The aberrations are characterised here by the amplitude of 
the Zernike polynomials representing the phase on each DIMM 
sub-aperture. They should not be confused with the aberrations 



Figure 9. Sequences of simulated spot pairs produced by a turbulent 
layer at 10 km with ro = 0. 1 m in a DIMM with perfect optics (top) and 
with a defocus of 1 rad (bottom). The same phase screen is used in both 
cases. 

of the feeding telescope. For example, a coma or spherical aber- 
ration in the telescope causes mostly astigmatism in the DIMM 
sub-apertures. Generally, the aberrations of both sub-apertures 
are not equal, but, to simplify, we assume here their equality. 

4.3 Propagation effects in DIMM 

The standard DIMM theory ('Martini 1 19871 : ISarazin & Roddied 
[1990) considers only near-field turbulence and neglects the prop- 
agation effects. DIMM is afected by three different phenomena 
related to propagation. 

Diffraction: part of the small-scale phase distortions are 
converted into amplitude fluctuations (scintillation) according 
to Eq. [T]and, as a result, the small-signal response of a DIMM 
slightly decreases with the propaga tion distance z. Th e decrease 
becomes noticeable for \/\z ^ D jTokoviniij|2002iJ) . as can be 
seen in Figs. [ToUl2l and[T4l 

Saturation. The second, even stronger effect is caused by 
the departure from the weak-perturbation theory and becomes 
important as soon as the scintillation index exceeds 0.1. Sim- 
ulations show that the response of a perfect DIMM to high- 
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Figure 10. Effect of propagation in a DIMM with D = 0.1m, iJ = 
0.24 m, and perfect optics. Monte-Carlo simulations with A = 0.5 ^tm 
and a single turbulent layer with different seeing located at the ground, 
at 5 km, and at 10 km. 
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Figure 11. Bias in the seeing estimated from the average longitudinal 
and transverse image motion (Eg. 112) as a function of the centroid win- 
dow radius S for near field (top) and 10 km propagation (bottom). The 
r.m.s. amplitude of the aberrations (coma ag, defocus 04 and astigma- 
tism 05) is 1 radian. Parameters: D = 0.1m, B = 0.25 m. 



altitude turbulence is non-linear, i.e. the coefficients Ki and Kt 
depend on the seeing (Fig.llOt. At 0'.'2 seeing, in the small-signal 
regime, the response is reduced only by diffraction. The satura- 
tion causes additional loss of response, depending on the seeing 
and layer altitude. On the other hand, the response to the near- 
ground turbulence remains constant even for a 3" seeing, despite 
strong distortions of the spots which split into several speckles. 
In practise, situations with > 0.2 are not uncommon and a 
DIMM is expected to "under-shoot". We could model this effect 
in the same way as we did for the MASS. In order to correct for 
the under-shoots, we need to know the turbulence profile. 

Aberrations together with propagation cause a complex 
bias considered in the next Section. It would be premature to 
correct DIMM for saturation and diffraction before the aberra- 
tion bias is addressed. 



4.4 Aberrated DIMM and propagation 

Small-signal response coefficients of an ideal and aberrated 
DIMMs are calculated with iA9\ and translated into the bias in 
the seeing e/eo derived from average longitudinal and transverse 
image motion, 



3/5 



(12) 



where the coefficients Kifl and Kt.o of an ideal DIMM are cal- 
culated by Eq. [TO] This bias is plotted in Fig. [TT] as a function 
of the centroid window radius 5 for two cases, turbulence at the 
ground (near field) and at 10 km distance. Both DIMM apertures 
have the same aberration. In the near-field case, the response is 
relatively stable against aberrations, as long as the centroid win- 
dow is wide enough, S > 2.5. 

When the spots are aberrated, the variance of the differen- 
tial image motion produced by a high turbulent layer depends on 
the type of the aberration, its amplitude, propagation distance, 
and the radius of the centroid window 5 (Fig. II It. 

The dependence of the DIMM bias on the defocus ampli- 
tude is further explored in Fig.[T2] Even for zero defocus, there is 
a small difference in the response between z — and z — 10 km 
because of the diffraction. The difference becomes larger for de- 
focused spots. The influence of a small defocus depends on its 
sign: on one side of the focus, the response decreases and DIMM 
"under-shoots", on the other side it increases. A defocus larger 
than 0.8 rad always causes a positive bias. These analytical cal- 
culations are confirmed by simulations. When ro = 0.5 m, the 
simulated response closely follows the small-signal curves. For 
To = 0.1m, the near-field response still matches the theory, 
but the 10-km response is reduced additionally by the saturation 
(Fig.El. 

A negative bias caused by the combination of propagation 
and small defocus may appear counter-intuitive: we expect that 
scintillation adds something to the image motion produced by 
the phase. Figure [T3] illustrates this effect from the geometric- 
optics perspective. Suppose that longitudinal image motion in a 
DIMM is created by a sinusoidal wavefront with a period ~ 2B. 
After propagation, the phase is reduced only slightly, but some 
amplitude fluctuations are created, positive where the rays con- 
verge and negative where they diverge (amplitude fluctuations 
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Figure 12. Bias of the seeing derived from the longitudinal (top) and 
transverse (bottom) image motion in a DIMM with D = 0.1m, = 
0.25 m, (5 = 5 for turbulent layers at the ground (dashed lines) and 
at 10 km (full lines), as a function of the Zemike defocus coefficient 
a4. The curves show analytical calculations for small signal, the results 
of simulations for tq = 0.1m and ro = 0.5 m are over-plotted with 
symbols. The dotted lines indicate Strehl ratio. 
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Figure 13. Illustration of the propagation effects in a DIMM (see text). 
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Figure 14. The dependence of longitudinal small-signal response coef- 
ficient Ki of a DIMM with D = 0.1m, B = 0.25 m, 5 = 5 on the 
propagation distance z for four values of defocus indicated on the plot. 
The behaviour of the transverse coefficient is similar. 



are proportional to the wave-front curvature). The gradient of 
the amplitude at DIMM apertures is correlated with the phase 
gradient (i.e. image motion). A perfect DIMM is insensitive to 
the amplitude fluctuations, but in a defocused DIMM the cen- 
troids are affected by the amplitude gradients. Depending on the 
sign of the defocus, the amplitude gradient either increases or 
decreases the measured image motion and thus creates a posi- 
tive or negative bias. 

A positive bias caused by strong defocus (and other aber- 
rations) has been confirmed experimentally by dir ect compar- 
isons between well-aligned and aberrated DIMMs jWang et"ail 
l2006h . The negative bias found here has not been suspected be- 
fore. A defocus of 0.3 rad rms (24 nm for A = 0.5 /^m) cor- 
responds to a Strehl ratio of 0.91, i.e. practically diffraction- 
limited spots. Yet, such a DIMM can have response coefficients 
reduced or increased by ~20%, and will produce seeing data bi- 
ased by ~ ±12% when most of the turbulence is concentrated 
at 10 km. Optical aberrations other than defocus are expected to 
cause a similarly complex bias. We conclude that the response 
of a DIMM to high-altitude turbulence is intrinsically inaccu- 
rate. The saturation makes things even worse, introducing de- 
pendence of the response on turbulence profile. 

We simulated the centroid calculation by thresholding nu- 
merically and found that its behaviour is very similar to the win- 
dowing method considered above. For turbulence at the ground, 
the response is insensitive to aberration, but slightly depends on 
the threshold (a bias of +10% in Ki and Kt for a threshold of 
0.1). The asymmetric dependence of response to high-altitude 
turbulence on the defocus (Fig.ll2t is also preserved. 

Figure[T4] shows how the small-signal response depends on 
the propagation distance z. For un-aberrated spots, the response 
slowly decreases with z. On the other hand, defocused spots 
show both positive and negative bias. A real DIMM instrument 
can either "over-shoot" or "under-shoot", depending on the state 
of its optics and the turbulence profile. 
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It is clear that the quality of the DIMM optics and focus has 
to be strictly controlled in order to get unbiased re sults. Mea- 
surement of the spot Strehl ratio (SR) suggested in (' Tokovininl 
|200 2a ) is helpful, with SR > 0.5 usually indicating acceptable 
quality dWang et alJlioO^ . However, the SR is reduced under 
poor seeing even in a perfect DIMM, so if the data with poor 
SRs are rejected, the seeing statistics will be biased. Yet an- 
other way to quantify the aberrations in a DIMM will be to 
take long exposures of defoc used spots and to analy s e them 
with the method described by iTokovinin & Heathcotel ilOOdl) . 
Needless to say that a DIMM working always out-of-focus (e.g. 
iBallv et al.|[r996.) will never measure the seeing accurately. 



4.5 Exposure-time bias 

Finite exposure time in a DIMM reduces the differential image 
motion, biasing the measured seeing to smaller values. The ef- 
fect can be quite stron g, as noted by Martin ( 1987). This bias was 
modelled in detail in ( Tokovinirill2002a() . Short (e.g. 5 ms) expo- 
sures reduce the bias, but do not eliminate it completely. If the 
data is acquired continuously, an extrapolation to zero-exposure 
can be done from the variances of time-binned signal, either by 
fitting an exponen tial curve to the d ependence of the variance 
on exposure time dWang et"aill2006l) or by using a simple lin- 
ear formula as in MASS or GSM dZiad et alJizOOOh . The linear 
extrapolation appears too strong, though. 

If the image sequence is not continuous, a method of inter- 
laced single and double exposures should be used. In this case, a 
"modi fied exponential correction" was developed in dTokovinii] 
If £i and £2 are the seeing values calculated with nom- 
inal and double exposure time, the de-biased seeing £o is esti- 
mated from 



eo ~ 0.5(ci£i + Ci''^e2), (13) 

where ci — (ei/ei)'^^'''- To reduce the statistical noise, the fac- 
tor ci is averaged (smoothed) over time and its average value is 
then used in il3\ to correct individual measurements. 



4.6 Centroid noise 

Even in the absence of atmospheric image motion, the measured 
centroids fluctuate because of the errors caused by the photon 
noise and detector readout noise. The errors of intensities lij 
are independent in each pixel and equal to the sum of readout 
and Poisson noise, expressed in the signal counts (ADU): 



(14) 



were R is readout noise in ADU, and G — the CCD camera con- 
version factor (gain) in e^/ADU. The influence of these errors 
on the calculated centroid is easily evaluated by differentiating 
Eq. [TT] and using independence of noise in each pixel. As the 



weights Wij take values of either 1 or 0, we simply restrict the 
summation to pixels where Wi,j — 1. The result is 

2 1 V^/ x2 2 



4 V(x,,.-c.)'(i?'+/.../G) 



(15) 



Formula 13 in jTokoviniij[2002al) contains a typographic error, 
reeled here 



Here Jo is the sum of intensities over pixels in the centroid win- 
dow. 

The sum entering in dlSt can be computed in advance if the 
centroid window has a well-defined size and the image profile is 
known. This is not the case when a thresholding method is used. 
However, even with thresholding the centroid noise of each indi- 
vidual spot can be evaluated with l |I5t during centroid computa- 
tion. Variations of the flux caused by scintillation or clouds can 
be accounted for as well. 

Obviously, the noise variance of both centroids in the 
DIMM has to be evaluated and subtracted from the measured 
differential variance before calculating the seeing with It 
would be wrong to express the noise in arc-seconds and sub- 
tract it later from the measured seeing, because the effects in 
the longitudinal and transverse directions are not the same. The 
noise depends on the stellar flux, hence subtracting a fixed "in- 
strumental noise" is not correct. Typically, the noise is small and 
its subtraction or not does not matter. However, the noise can 
significantly bias DIMM results under good seeing or for faint 
stars. The flux in the interlaced single and double exposures is 
different, hence if the noise variance is not properly subtracted, 
it will bias (over-estimate) the correction l |13l l. 

To reduce the noise, especially the term related to R, the 
smallest number of pixels must be used, favouring a narrow cen- 
troid window with 5 ~ 1 or, equivalently, a high threshold. How- 
ever, we have seen in Sect. 14.21 that such choice leads to a very 
uncertain response coefficient in presence of even small aber- 
rations. The noise is usually much less of a problem than the 
aberrations, so a wide centroid window with S ^ 2.5 or a low 
thres hold are prefera ble, contrary to the recommendation given 
in aTokoviniij|2002ah . 

Both the measurement noise and the noise caused by the 
scintillation are isotropic and affect longitudinal and transverse 
image motion equally. The effect on the calculated seeing, how- 
ever, will be different because Ki > Kt- If the DIMM signal 
contains a significant noise component, the "transverse" see- 
ing will always be larger than the "longitudinal" seeing. This 
effect could be mistakenly interpreted as a manifestation of a 
non-Kolmogorov turbulence spectrum. 



5 CONCLUSIONS 

This work complements previous studies of two methods to 
measure seeing - DIMM and MASS - and focuses on the 
achievable accuracy. The "seeing" itself cannot be defined very 
accurately, being a model-dependent parameter of a random and 
non- stationary process. Taking aside this caveat, we investigate 
potential instrumental biases by simulating both turbulence and 
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the instruments numerically. The true seeing is known exactly, 
permitting to calibrate the methods on the absolute scale. 

Our conclusions and recommendations can be summarised 
as follows: 

• The departure from the weak-perturbation theory affects 
the MASS method seriously, but it is correctable for a not-too- 
strong scintillation, s\ < 0.7. 

• The effects of finite exposure time and inner turbulence 
scale in MASS can be neglected. 

• The DIMM method is very robust and tolerant to aberra- 
tions for the near-ground turbulence, provided that the centroid 
calculation uses a large enough radius, > 2.5X/D. 

• Optical propagation causes a non-linear response of DIMM 
(saturation), not present in the near-field case. 

• Combination of propagation and aberrations in a DIMM 
leads to a complex bias. By controlling the DIMM optical qual- 
ity (Strehl ratio > 0.5), we can keep the seeing bias to within 
±12%, but making the control tighter appears impractical. 

• Centroid noise in a DIMM should be computed and sub- 
tracted from the measured variance. To do so, the detector read- 
out noise and gain must be known. 

• The exposure-time bias in a DIMM should be corrected us- 
ing one of several known recipes. 

The sensitivity of the DIMM to propagation and aberra- 
tions comes as a surprise, although neither of these effects was 
accounted for by the standard, near-field DIMM theory. The 
bias on high-altitude seeing is so complex that removing it com- 
pletely seems unrealistic. 

It has been demonstrated that two identical DIMM instru- 
ments with good optics can give seeing meas urements concor- 
dant to within few percent jWang et"alll2006l) . Are these mea- 
surements accurate to the same level? Not necessarily. A seeing 
of 1" coming from 10 km could be measured by both instru- 
ments with a bias of 0.83 (Fig. |10| (. whereas the same seeing 
originating near the ground will be measured correctly. 

The response of a DIMM depends on both instrumental fac- 
tors and observing conditions. Two different DIMMs can agree 
on one night (e.g. seeing dominated by low layers) and dis- 
agree on another night or at another site. Two identical DIMMs 
can be biased in a different way at two different sites. Thus, 
inter-comparison between DIMM (or MASS) instruments can- 
not guarantee that they are accurate. Their mutual agreement 
is a necessary, but not sufficient condition. Only a careful con- 
trol of biases can ensure accurate seeing data. However, instru- 
ment inter-comparisons are useful for debugging and checks and 
should be pursued whenever possible. 

A combination of MASS and DIMM in one instrument has 
stimulated this research. These instruments, when properly cal- 
ibrated, agree very w ell for a seeing dominated by high layers 
jKornilovet"ai]|2007l) . This is a triumph of the optical propaga- 
tion theory enabling us to interpret both scintillation and image 
motion with a common model and a single parameter, tq. At the 
same time, the agreement between two instruments based on dif- 
ferent principles and with different biases is a strong argument 
that both are accurate, i.e. measure the seeing on the absolute 
scale. 
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APPENDIX A: RESPONSE OF A CENTROID 
ESTIMATOR IN THE PUPIL PLANE 

The centroid signal c for each DIMM spot is obtained from the 
weighted PSF I{a) as 



da I{a) M{a) 



(Al) 



where a is a 2-dimensional vector of angular coordinates, M{a) 
is some function, mask, and Jo is the total intensity (flux). All 
integrals are in infinite limits and exist because all functions 
are supposed to have limited support. The PSF is not neces- 
sarily an ideal one, but may include some aberrations. For- 
mula JAU is rather general and applies to many situations, e.g. 
curvature sensing. For centroid calculation, we need a mask 
M{a) = axw{a) to match Eq.[TT] We keep iAU in a general 
form useful for other applications. 

Let X be the coordinate vector in the pupil plane. The com- 
plex amplitude of the initial un-perturbed field at the pupil is 
U{x). It includes the pupil function (possibly with aberrations) 
and is normalised arbitrarily. Suppose that the amplitude U is 
changed by a small phase aberration ip{'K) and a small log- 
amplitude perturbation x{^) ^rid becomes U{x) 6"^^''''+'^''''' . 
What would be the change of the signal Ac caused by this aber- 
ration? 

We find a small change in the signal b y linearising the 
known expression of the OTF dOoodmarJl 19851) 



/(/) = Jo"' J d^x U{x)U*{x + Xf) 



(A2) 



with respect to small p erturbations y <C 1 and x ^ 1> so-called 
PSF Taylor expansion jPerrin et al.ll2003h : 

A/(/) = lo' J d^xU{x)U*{x + Xf) (A3) 
X [Mx)+x{x)-if{x + \f)+x{x + Xf)] 
The signal increment is 



Ac= / dV A/(/)A/(/). 



(A4) 



We put l |A3t into l lA4t and re-group the terms. The first term in 
the square brackets containing iipix) leads to 

Aci = i(A^7o)"^ j d^x j d^x' U{x)U* {x + x') 

X M{x /X) ip{x), (AS) 

where x' — Xf. The 3-rd term containing —i(fi{x + Xf) leads to 
the complex-conjugate of the expression l |A5t with inverse sign. 
Collecting all 4 terms, we write the result as 

Ac = y" d^x F^{x) ifi{x) + J d^x F^{x) x{x), (A6) 

where the filter functions F are 

F^{x) = Im[A{x)], F^{x) = Re[A{x)] (A7) 



and 



A{x) = 2 (A^/o)"^ U{x) J d^x' Vix + x') AI{x'/X) (A8) 

The response is independent of the normalisation of the ampli- 
tude U because it is divided by the flux. 

This result contains an implicit assumption that the fluctua- 
tions of the denominator 7o can be neglected. This is not always 
true. The signal has the form c — A/B, hence its fluctuations 
are Ac = A A/B — c{AB/B). The fluctuations of the denom- 
inator AB can be neglected if the average signal c — (they 
will be a second-order term then). This condition is enforced 
by the choice of Af(a) — ax'w{a) appropriate for a windowed 
centroid. 

If the mask AI is correctly dimensioned to compute the cen- 
troid c in pixels and the angular size of the pixel is p, then it 
follows that the longitudinal response coefficient Ki is 



dV 



Ki = {pD/Xf {ro/Df^'' 

X <^>^(/)[2sin(^B^)]^ 



(A9) 



where both centroid measurement direction and baseline are par- 
allel to the X axis. To calculate the transverse response Kt, we 
replace fx with fy. To take into account both propagation and 
sensitivity to scintillation, we make a replacement 



\FAf)\' ^ \FM)cos{nXz\f\^)+FM)sin{^^^\fi 



.(AlO) 



This modification automatically accounts for the correlation be- 
tween phase and amplitude. 

This paper has been typeset from a T^i/ KI^ file prepared by 
the author. 
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